AFOSR^7-l$*7 

CALIFORNIA  INSTITUTE  OF  TECHNOLOGY 

Pasadena,  California 

ANTENNA  LABORATORY 

Regoyt  No.4i 


Electronics  Division  f'This  c'oainont  has  Eccn  e-rrovi; Y\ 

.  j  lor  taiblv  tnv^  re1-;  I  ;  ] 

Air  Force  Office  of  Scientific  Research  |  d^^ibaiicm  k.  nabuLtci,  \ 

Contract  No.  A F4 9 ( 638^) -1 266  ~ 

August  1967 


C'AMFOPSIA  :*iSTIT',T£  OF  TECFbblVJY 
Pasadena,  California 


Antenna  Laooratsr;' 
TECHNICAL  REPORT  UO.  1*1 


POTENTIAL  OK  a  CHARGED  CY1IV1'-5 


by 


Jcnn  bar. 


A FOSE  67-1667 


This  research  was  supported  by  the 
U.  S.  Air  Force  tnrougn  tr.e  air 
Force  Office  of  Scientific-  Research 
of  the  Office  of  Aerospace  Research, 
unaer  Contract  AFb9\b38)~i266. 
Reproduction  in  whole  or  in  part  is 
permitted  for  any  purpose  of  the 
United  States  Government-. 


August  1967 


2. 

ABSTRACT 

A  method,  based  on  complex  variable  theory,  is  presented 
Tot  solving  potential  problems  involving  conducting  solids  of 
revolution.  Attention  is  paid  excliisively  tc  the  right  circular 
cylinder,  both  because  this  problem  possesses  a  high  degree  of 
symmetry  which  brings  about  certain  simplifications ,  and  because 

\ 

i  comparison  can  be  made  with  previous  results.  The  generalisation 

e 

from  the  cylinder  to  other  axially  symmetric  geometries  should  . 

|  become  obvious.  This  method  depends  on  the  establishment  of  an 

integral  representation  for  the  potential  function,  which  leads 
\  .  to  the  formulation  of  the  problem  in  terms  of  a  pair  of  coupled 

Fredholm  integral  equations  of  the  first  1'ind.  As  an  illustrative 
example  these  equations  are  solved  numerically  in  the  lowest 
approximation  and  the  capacitance  of  the  cylinder  is  calculated. 
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1.  Introduction 

I r  '’-his  work  ve  present  a  method  for  solving  the  electrostatic 

problem  of  a  charged*  right  circular,  conducting  cylinder.  This 

1  2 

problem  has  beer  treated  previously  by  faythe  *  ,  using  a  method 
which  depends  cn  tl;c  following  very  important  property  of  an  axially 
symmetric  potential  function*  namely*  the  value  of  the  potential  every¬ 
where  is  known  once  its  value  on  the  synmetry  axis  is  determined.  Our 
method  also  depends  on  this' property *  but  in  quite  a  different  way. 

Here  wc  use  it  t o  “est=bi5  sh  an  Integra],  representation  for  the  poten¬ 
tial  functicn  which  has  a  very  simple  form.  The  problem  is  then 
formulated  in  terms  of  a  pair  of  integral  equations  which  are  well 
suited  for  numerical  integration.  Inasmuch  as  Snqrthe  has  published 
numerical  results  of  high  accuracy,  it  is  not  our  primary  intention  to 
reproduce  these  solutions  in  great  detail.  Ho  heavy  computation  is 
therefore  attempted.  Our  interest  is  mainly  methodological.  It  is 
hoped  that  this  method  can  be  applied  with  equal  ease  to  other  solids 
of  revolution. 

2.  Whittaker's  Solution 

Denote  the  radius  of  the  cylinder  by  a  and  its  length  2b  . 

Let  us  set  up  a  cylindrical  coordinate  system  such  that  the  origin  and 
the  z-axis  coincide  with  the  center  and  the  synmetry  axis  of  the  cylin¬ 
der.  Thus  the  cylinder  is  defined  by  ths  intersecting  surfaces  o  -  a 
and  z  =  +b  .  The  potential  function  we  seek  is  a  solution  of  the 
axially  symmetric  Laplace  equation 


-  , 
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This  is  a  special  case  of  Uhit.taker's  solution  cf  Laplace's  equation  . 
It  gives  a  prescription  for  obtaining  the  off-axis  value  of  the  poten¬ 
tial  in  terns  of  the  value  on  the  upper  axis  provided,  of  course,  that 
the  latter  is  as  analytic  function.  Ve  note  that  for  -b  <  z  <  b  , 
V(0,z)  is  undefined.  Its  value  there  must  be  obtained  by  analytic 
continuation. 


The  potential  function  Y(p,z)  has  another  integral  representation 
in  the  form  of  the  particular  solution  of  Poisson's  equation: 


V(p,z)  * 


(3.1) 


where  o  is  the  surface  charge  density,  R  the  distance  from  a  point  on 
the  cylinder  to  an  observation  point,  and  the  integral  is  taken  oyer  the 
surface  of  the  cylinder.  For  points  on  the  axis 


(3.2) 


The  three  integrals  represent  contributions  from  charges  on  the  top,  the 
bottom  and  the  side  of  the  cylinder  respectively.  From  symmetry  ve  have 


o(p',b)  *  o(p',-b) 

o(a,z')  *  o(a,-z') 


(3.3) 


The  o's  are  continuous  and  bounded  except  at  the  edges  p  =  a  , 

2|  *■* "  — _ 

*  *  ^  §  where,  according  to  the  edge  condition  %  they  diverge  only 
-1/3 

like  6  ,6  being  the  distance  Tree  the  edge.  Thus  (3.2)  defines 

V(0,z)  as  an  analytic  function  of  z  .  Clearly  V(0,z)  is  real  for  z 
real  and  tends  to  zero  like  z_1  as  |z|  ■*  •  . 

The  singularities  of  V(0,z.)  can  be  read  off  from  (3-2)  r  The  first 
integral  has  a  branch  cut  from  b-ia  to  b+ia  ;  the  second:  integral  hoe 
a  cut  from  -b-ia  to  -b+ia  ;  and  the  third  integral  has  two  cuts  from 
-b+ia  to  b+ia  respectively.  These  branch  cuts  of  V(0,z)  from  a  rec¬ 
tangle  in  the  z-plane  which  is  identical  Xo  the  section  of  the  cylinder 
by  a  plane  containing  the  symmetry  axis.  It  can  be  shown  that  this  prop¬ 
erty  of  V(0,z)  is  true  of  all  solids  of  revolution. 

The  symmetry  of  the  problem  in  the  physical  space  implies  a  certain 
symmetry  in  the  z-plane.  We  readily  deduce  from  (3.2)  and  ( 3* 3)  that  the 
absolute  values  of  V(0,z)  and  V-(0,-z)  are  equal.  If  ve  analytically 
continue  V(0,z)  from  the  point  z  to  the  point  -z  ,  we  must  choose  a 
path  which  avoids  the  branch  cuts  near  the  origin.  Then  during  the 
process  all  the  square  roots  in  (3.2)  change  sigh.  We  conclude  that 
whereas  the  value  of  the  potential  on  the  symmetry  axis  is  an  even  function 
of  z  in  the  physical  space,  it  is  an  odd  function  in  the  complex  z-plane: 

V(0,-z)  =  -  V(O.z)  (3»!*) 

In  general,  the  value  of  V(0,z)  or*  the  positive  real  axis  of  the  z-plane 
coincides  with  the  physical  value  of  the  potential  on  the  upper  symmetry 
axis,  while  the  value  on  the  negative  real  axis  differs  in  sign  from  the 
physical  value  on  the  lower  symmetry  axis.  This  is  easily  seen  from  (2.2) 


?. 
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Equating  the  real  parts  of  (3*5)  and  (3-T}»  ve  get 


Vffjr)  -  R«  -t  \ 

VI  i 


o(o'.b)  c '  do1 


o  yp'2  -  / 


a 

-*  i 


o(p'  ,b)  p'  dp* 


(3.8) 


d  ,  H2  T 
|yi  vp  -  y 


From  this  result  ve  see  that  f(y)  is  continuous  in  the  interval 
-a  *  y  -  a  . 

A  similar  integral  representation  for  the  second  integral  in  (3.2) 
is  obtained  by  putting  -b  for  b  in  (3.5) : 


a 

ir  \ 

n  * 


o(p’.*-o-)  p'  dp' 


vQ  r  f(y)  <*y 

¥  i  -b  ♦  iy  - 
-a  * 


(3.9) 


0  Jp'2*lz*b)2 

The  third  integral  in  (3.2)  has  two  branch  cuts  and  is  represented 
as  follows: 


1 

2e 


o(a,z' )  a  dz' 


b+ia  &,U)  dC  b"fia  &2U)  d£ 


— V)  Jar  -  (z  -  z*)‘ 


[  glU;  lo  f 

"  "  ’  J  5  -  z  "  *  J 


-b+ia 


-b-ia 


C  -  z 
(3.10) 


where  g^(0»  g9(£)  are  real  density  functions.  The  left-hand  side  is 
real  for  z  real;  hence 


S1(C)  =  ggU)  =  g(S) 


(3.11) 


and  we  have 
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P 

±  1 
n  * 


o(a,s‘)  adz' 


b  i  J2  >  _7v 

-d  yfe  T  i*-‘  ; 


v&-j  g(x)  dx  VQ  |  g(x)  dx 
*  lx*  ia  -  z  *  !  x  -la  >  2 


(3.12) 


By  (3.4),  g(x)  is  an  even  function 


g(-x)  *  g(x) 


(3.13) 


It  can  also  be  shown  that  g(x)  is  continuous  in  the  interval  -b  *  x  ^  b 
Collecting  all  the  results  ve  have  the  following  integral  representa¬ 
tion  for  V(0,z) ;  - 


V  f  f(y)  4*  V  >  f(y)  ay 

v(o,i) - 1  -  — >■ .  -f  !  — — — 

1  b  ♦  iy  -  z  £  -b  f  iy  -  z 


Is.  }  8(>)  i  }  g(x)  ** 


i.  f 

W  i  x  + 


I3.1V) 


ia  -  z 


4  x-ia  -  z 


where  f(y)  and  g(x)  t>re  a  pair  of  real,  continuous,  even  functions. 


4.  An  Integral  Renrasgntat ion  for  V'o  ,z) 


We  replace  z  by  z  *  ip  cos  a  in  (3.14),  substitute  the  result 
in  (2.6),  and  interchange  the  order  of  integration.  The  integrals  over 
a  have  the  form 


1! 

1  f  — 

’  l  A* 


(4.1) 


B  cos  a 


where  A  and  B  are  cou^iex  numbers .  By  an  identity  of  Jacobi  (4.1)  is 


equal  to 


1 


(4.2) 
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where  "the  sign  of  the  square  root  is  chosen  to  satisfy  the  inequality 


|A  -  Va2-  B2|  <  |  B  | 


»e  therefore  obtain  the  integral  representation 
V  ?  f(y)  dy 


V(p 


-a  ^z-b  -  iy)2  ♦  p2 


Vo  |  f(y )  dy 


-a 


V  uf 
* 


g(x)  dx 


♦  i  f 


|/(z  +  b-i  r>2+  p2 


g(x)  dx 


-b  ^z-x-ia}2+  p2  -b  |j( z-x+ia )2+  p2 


(fc-3) 


(fc.M 


For  p,z  large  ana  positive,  the  signs  of  tne  square  roots  are  chosen  to 
■ake  the  real  parts  positive.  For  other  ranges  of  value  the  signs  ere 
deterained  by  continuation.  It  can  be  shown  that  the  real  parts  are  taken 
to  be  positive  throughout .  Using  the  evenness  and  reality  of  f(y)  and 
g(x),  we  can  also  write  (h.k)  in  the  fora 


V(p,z)  = 


2V  7 

"*“1 


f(y)  dy 


^z  -b  -  iy)2+  p2 

g(x)  dx 


2V  7 


2V  7 


-b  |/(z  -  x  -  ia)2+  p2 


f(y)  ay 


0  j^z  +b  -iy)2+  p2 


(«*.5) 


The  capacitance  is  given  by 


16Sq  [  f(y)  dy  +  |  g(x)  dx j 


(U.6) 


The  boundary  conditions  for  V(p,z)  are 


V(p,+b)  =  V 


0  <  p  <  a 


(**.7) 
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V'(ft,z)  =  V  ,  -b  <  z  <  b 

c 

Combining  {1* - 5 ) »  (4--7J  and  (4.8)  ve  have 

f(y  ;  dy 

/(2b  -  iy)2+  b2 


1  [ 

*  !  /  x - ~  ~  1 


n  /~2”  2 

0  1/0  -  y 


(4.8) 


+  -  Re 
ir 


g(x)  dx 


-b 


y<x  -b  +  ia) 


-  1  , 


0  <  p  < 


(4.2) 


2a  2 

♦  o 


f  f(y)  dy  of  f(y)  dy 

V  Re  j  - - +  f  Re  j  . . . 

G  J( b-z  +  iy)2+  a2  0  J(b  +  z-  iy)2+  a2 

2  r  g(x)  dx 

+  £  Re  j  -  =  1  ,  -b  <  z  <  b  (  4 . 1C ) 

11  ■*  r~  2  2 

-b  J[x -  z  +  ia)  +  a 

These  are  a  pair  of  coupled  Fredholm  integral  equations  of  the  first  kind 
which  summarize  the  formulation  of  the  problem.  We  note  that  in  the 
usual  integral  equation  formulation  of  the  problem  we  start  with  the  inte¬ 
gral  representation  (3-1)  and  impose  the  boundary  conditions  (*4.7)  and 
(4.8).  In  this  way  we  also  obtain  a  pair  of  coupled  integral  equations— in 
o(p',b)  and  o(a,z')  .  But  here  the  kernels  are  elliptic  integrals. 
Moreover  the  surface  charge  densities  are  known  to  diverge  at  the  edges. 

5.  Two  Special  Cases 

Before  proceeding  to  discuss  the  solution  of  (4.9)  and  (4.10)  we 
shall  first  examine  a  few  special  cases. 
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Case  I:  A  Circular  Disk 

Letting  b  =  0  in  (n.9)  we  obtain  the  limiting  case  of  a  charged 
circular  disk: 

4  (p  f(y)  dy 

. . .  1  .  0  <  p  <  a  (5.1) 

2  2 

-  y 


This  is  a  special  case  of  the  Abel  integral  equation 
jP  f(y)  dy 
0 


J7: 


-  F(p) 


0  <  p  <  a 


(5.2) 


The  solution  is 


f(y)  =  f 


-  f 

dy  J 


*r  F(p)  p  dp 


ft- 


0  <  y  <  a 


(5-3) 


Therefore  we  have 


and 


V(p,z)  = 


f(y)  =  | 


dy 


-a 


J{ z  -  iy)2  +  p* 


(5.«0 


V  z  -  la  + 

o  .  . 

—  i  £n  - 


y/(z  -  ia)2+  p2 


r  22 

z  +  ia  +  ^(z  +  ia)  +  p 


2V 


°  ,  -1 
tan 


2  a 


-|l/2 


|  2^  2  2  [72 “  2  2.2  ,22 

Lp  +  z-  a+  |/(p  ♦  z  -  a  )  +  4a  q.  J 


(5-5) 
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Our  method  of  solution  of  to  is  prcblea  is  closely  related  to  that  of 
Heir.s  and  XacCaay'*. 

Case  II:  Two  Coaxial  Circular  Disks 


Letting  g(x)  »  0  in  (t.G)  we  have 


rp  f(y)  dy 


k  I?  2" 

0  i/p  -  y 


-I  f 

*  J 

-a 


f(y)  dy 


(2b-  iy)2+  c2 


0  <  s  <  a 


(5.6) 


where  the  upper  and  lower  signs  correspond  to  the  cases  of  equally  and 
oppositely  charged  disks  respectively.  If  we  treat  the  right-hand  side 
of  (5-6)  formally  as  a  known  function  of  c  ,  this  is  an  Abel  integral 
equation.  Upon  applying  (5-3)  we  get 


*»>  ±f  I 


f(y’)  dy’ 


(5.7) 


JO  o 

a  +  (y’~  y) 


which  is  a  Fredholm  integral  equation  of  the  second  kind.  This  equation 
was  first  derived  by  Love^ ,  who  also  obtained,  by  a  different  method,  an 
integral  representation  similar  to  (k.U). 


numerical  Solution  of  the  Integral  Equations 


We  are  unable  to  solve  the  equations  ( - 9 )  and  (Ij.10)  exactly. 
This  should  become  quite  clear  in  the  last  section  where  we  could  not 
even  obtain  an  exact  analytic  solution  for  the  special  case  of  two 
coaxial  disks.  However,  the  form  of  these  equations  is  simple  enough 
for  them  to  be  integrated  numerically.  A  numerical  solution  usually 
gives  f(y)  and  g(x)  as  polynomials  of  y  and  x  respectively. 


Ik. 


Alternatively  ve  can  approximate  f(y )  and  g(x)  by  polynomials 
from  the  start: 


w  -  !  vf’2" 

n=0 

«<*>  -  I  (6.D 

*■0 


Substituting  these  expressions  into  (4.9)  and  (4.10)  and  carrying  out 
the  integration? .  ve  obtain  two  functional  equations  involving  the 

3  ♦  M  +  2  coefficients  A.,  and  .  In  these  equations  ve  set  0 

2m 

successively  equal  to  3+1  values  in  the  range  0  <  0  <  a  and  z 
equal  to  M+l  values  in  the  range  0  <  z  <  b  •  In  this  way  ve  obtain 
3  ♦  M  ♦  2  linear  algebraic  equations  for  the  determination  of  the 
unknovn  coefficients. 

As  an  illustration  ve  use  the  lovest  approximation  by  taking  only 
the  first  terms  in  (6.1): 


i’(y)  s  Aq  ,  g(x)  =  Bq 


(6.2) 


Putting  0 


z 


respectively  in  (4<9)  and  (4.10)  ve  get 


a,  .A  +  al0B 
li  a  12  0 


a21ft0  *  S22Bo 


In  terms  of  the  ratio  A  =  b/a 


,  .  2  -1 
a.,  =  1  +  —  tan 

11  TT 


=  1 

=  1 

,  the  coefficients  are  given  by 

1+1 


(6.3) 


2A  +  a 


The  capacitance  is  given  by 

C  =  l6e  (A  +  AB  )  a  (6.6) 

o  o  o 

Letting  A  =  0,  1,  2,  U  successively,  ve  obtain  the  following 

results : 


A 

A 

o 

B 

o 

C  in  farads 

Sraythe's  value 

0 

.500 

__ 

.708  x  10-I°a 

.708  x  icfl0a 

i/i» 

.578 

.310 

.928  x  10-10a 

.922  x  10-I0a 

1/2 

.625 

.290 

1.0'-  x  10-I0a 

1.07  *  10-10a 

1 

.63*4 

.330 

1.36  x  10-10a 

1.33  x  10-10a 

2 

.588 

-3Ul 

1.79  *  10-10a 

1.75  x  I0~10a 

k 

.531 

.312 

2.52  x  10-I0a 

2.U7  x  I0”10a 
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In  this  approximation  the  capacitance  differs  from  Scythe ' s  value  by 
about  2%. 


However,  if  we  use  the  above  solution  to  compute  the  surface  charge 
density,  we  find  that  our  value  deviates  quite  seriously  from  the  result 
of  Snythe.  This  is  not  surprising  in  view  of  the  crudity  of  our  approxi¬ 
mation  and  the  fact  that  the  charge  density  is  a  local  quantity.  The 
capacitance,  on  the  other  hand,  is  a  global  quantity  and  is  not  very 
sensitive  to  the  number  of  terms  we  take  in  (6.1)  to  represent  f(y)  and 
g(x) .  Taking  move  terms  will  certainly  increase  the  accuracy.  Neverthe¬ 
less,  we  will  not  go  further  into  elaborate  calculations. 
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II  SUPPUMCN1ARV  NOTIJ 


TECH  CTHER 


•J  A0STHACT 


A  method,  based  on  complex  variable  theory,  is  presented  for  solving 
potential  problems  involving  conducting  solids  of  revolution.  Attention  is  paid 
exclusively  to  the  right  circular  cylinder,  both  because  this  problem  possesses 
a  high  degree  of  symmetry  which  brings  about  certain  simplifications,  and  because 
comparison  car;  be  made  with  previous  results.  The  generalization  from  the  cylinde 
to  other  axially  symmetric  geometries  should  become  obvious.  This  method  depends 
on  the  establishment  of  an  integral  representation  for  the  potential  function  which 
leads  to  the  formulation  of  the  problem  in  terms  of  a  pair  of  coupled  Fredholm 
integral  equations  of  the  first  kind.  As  an  illustrative  example,  these  equations 
are  solved  numerically  in  tne  lowest  approximation  and  the  capacitance  of  the 
cylinder  is  calculated. 
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